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Abstract. We present a theory of the dielectric response of a solution containing 
large solutes, of a nanometer size, in a molecular solvent. It combines the molecular 
dipole moment of the solute with the polarization of a large subensemblc of solvent 
molecules at the solute-solvent interface. The goal of the theory is two-fold: (i) to 
formulate the problem of the dielectric response avoiding the reliance on the cavity- 
field concepts of dielectric theories and (ii) to separate the non-additive polarization 
j\ of the interface, jointly produced by the external field of the laboratory experiment 

• i-h . and the solute, from specific solute-solvent interactions contributing to the dielectric 

^_j signal. The theory is applied to experimentally reported frequency-dependent dielectric 

^S^ • spectra of lysozyme in solution. The analysis of the data in the broad range of 

, £~ ( frequencies up to 700 GHz shows that the cavity field susceptibility, critical for the 

theory formulation, is consistent with the prediction of Maxwell's electrostatics in the 

frequency range of 10-200 GHz, but deviates from it outside this range. In particular, 

it becomes much smaller then the Maxwell result and shifts to negative values at 

f^ ■ small frequencies. The latter observation implies a dia-electric response, or negative 

^D , dielectrophoresis, of hydrated lysozyme. It also implies that the effective protein dipole 

recorded by dielectric spectroscopy is much smaller than the value calculated from 



protein's charge distribution. We suggest an empirical equation that describes both 
the increment of the static dielectric constat and the decrement of the Debye water 



peak with increasing protein concentration. It gives fair agreement with broad-band 
dispersion and loss spectra of protein solutions, but misses the (5-dispersion region. 
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1. Introduction 

Dielectric spectroscopy is a linear response technique, monitoring the dynamics of the 
dipole moment of a macroscopic sample of a polarizable material [TJ |2] . While it is highly 
sensitive and provides wealth of information about the dynamics of polarization modes 
active in a medium, the interpretation and the assignment of the observed relaxation 
processes often require theoretical approaches. 

The standard theoretical tool to study mixtures is the Maxwell- Wagner theory [21 E] 
and its modifications, also in terms of effective- medium approaches [3]. All these theories 
assume that macroscopic dielectric constants can be assigned to all components of the 
mixture. This often becomes a significant oversimplification when highly heterogeneous 
solutes of nanometer dimension, such as hydrated proteins, are involved |6]. The 
description of the polar response in terms of the molecular charge distribution is more 
accurate for these solutes [7J. Given the length-scale of the external field variation in 
a typical dielectric or light-absorption experiment, the overall charge and the dipole 
moment are the two main multipoles to consider (8]. 

Once a dipole is assigned to a protein, one might assume that standard models of 
dipolar liquids [TU], involving the statistics and dynamics of molecular dipoles, can 
be directly extended to study protein solutions. One has, however, to recognize that 
proteins, and other solutes of similar dimension, possess an extended interface with a 
molecular solvent, such as water, which is absent in the case of mixtures composed 
of molecules of comparable size. The interface of a hydrated protein involves a large 
number, ~ 300 — 500, water molecules only in the first hydration layer. Given that 
the perturbation of the water polarization propagates at least into the second hydration 
layer [HI [12] , the actual size of the protein- water interface is significantly larger. 

These new physical realities pose the requirement to develop new theoretical 
approaches to describe the polar response of protein solutions. The key question for 
this development is how to extend the classical theories of polar response of molecular 
dipoles into the realm of large solutes with an extended interfaces. The key parameter 
for the development of dielectric theories is the Onsager cavity (or directing) field [13] 
producing the torque on the molecular dipole when the macroscopic sample is placed in 
an external electric field [2]. The standard result of the classical theories is that the 
field of external charges E is screened by the polarization of the interface to the cavity 
field [H EH] 

E c = — - — E , (1) 

where e s is the dielectric constant of the solvent. This cavity field then directly leads 
to the Onsager mean-field [16] equation for the dielectric constant and, when mutual 
short-range correlations of dipoles are included, to the Onsager-Kirkwood relation [13] . 
The problem one faces in an attempt to describe a mixture of nanometer-size solutes 
with a molecular solvent is that there is no analog of either of these two equations. 
The fundamental line of inquiry here is whether one can extend equation ([T|), or its 
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analog, to such mixtures, or a new set of rules is required. This question is nontrivial to 
answer, even though some initial computer simulations indicate that, indeed, polarized 
nanoscale interfaces follow rules different from those established for cavities carved in 
dielectrics [T71 HH [12] . The simulations are however limited by the nanosecond range of 
time-scales. The question of what is the polar response of a nanoscale interface at low 
frequencies remains therefore open. 

This study aims to address this question by analyzing recent dielectric data obtained 
for solutions of lysozyme in water [191 EO] . We first develop a general formalism that 
does not anticipate any particular solution for the local field acting on the protein dipole. 
As a result of this analysis, we arrive at a surprising conclusion that the hydration layers 
of the protein screen its dipole even more substantially than anticipated by the standard 
result for a dielectric cavity given by equation (CQ). 

We start with introducing the polarization of the solute-solvent interface by the 
combined effect of the external electric field and the solute dipole moments. This 
interfacial polarization integrates into an interface dipole moment, which is assigned to 
each solute even in the absence of its own dipole. This development leads to the equation 
for the dielectric constant of an ideal solution of dielectric voids inside the polar liquid. 
We show that this equation is quite useful in describing the high-frequency dielectric 
response of a real solution, when the relaxation of the solute dipoles is dynamically 
frozen. We then proceed to a mixture of polar solutes in a polar liquid. Here, the 
cross-correlations of the solute and solvent dipoles [21] are expressed in terms of the 
cavity-field susceptibility, which can take different forms depending on the microscopic 
structure of the water layer interfacing the solute [18] . 

2. Dielectric Response of a Mixture 

Dissolving a polar solute in a polar solvent leads to two distinct effects on the response 
of the medium to an external electric field. The first effect is the exclusion of the 
solvent from the volume of the solute. The second effect is the response of the charge 
distribution within solute to the orienting torque of the external field. The two effects are 
entangled in the polarization of the interface by the solute charges and by the external 
field. However, their contributions to the overall dielectric response of the solution can 
be separated in the frequency domain. Since they also originate from distinct physical 
interactions, repulsive expulsion on the one hand and electrostatic interactions on the 
other, we start with considering the effect of the solute excluded volume and then add 
the contribution of the solute dipole moment to the dielectric response of the solution. 

2.1. Non-polar Solutes in a Polar Solvent 

Excluding the solvent from the solute volume creates the solute-solvent interface. 
From the standard viewpoint of dielectric theories, any interface carries an interfacial 
polarization when the solution is placed in a uniform field of external charges (capacitor 
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plates of the dielectric experiment). The polarization of the interface is described in 
the Maxwell theory of dielectrics by the surface charge density [T5]. It is given as the 
projection of the dipolar polarization of the dielectric P^ at the dividing surface S on the 
outward normal to the surface n^. The surface charge density then becomes op = n-Ps- 
This charge density integrates to a dipole moment of the interface 

M> nt = fr s a P (r s )dS, (2) 

Js 

where the surface integral is taken over the closed surface S enveloping the solute. 

The interface dipole polarizes the surrounding dielectric by its own electric field such 

that the inhomogeneous Maxwell field E(r) around the solute is a sum of the uniform 

Maxwell field of the external charges e ^Eq and the dipolar field of the polarized interface 

E(r) = e ; 1 E + ^T(r-r J )-M| ) nt . (3) 

5 
Here, T(r — r,) is the dipolar tensor describing the electric field at point r inside the 
solvent by a point dipole placed at ry; the sum runs over Nq solutes with coordinates 

The Maxwell field E(r) polarizes the liquid, with the resulting local inhomogeneous 
polarization P(r) = (47r)~ 1 (e s — l)E(r), decaying to the homogeneous polarization P of 
the external charges far from the solute-solvent interface. The overall dipole created in 
the solution is the integral of P(r) over the volume Q occupied by the solvent 

M mix = /p(r)dr. (4) 

Jn 

Here, subscript "mix" identifies the solvent-solute mixture. Assuming that the 
interfacial dipoles of solutes are independent of each other, one gets [32] 

M rnix = M hom - N Q Q P - (2/3) (e, - 1)N M™\ (5) 

Here, Mh om — VP is the dipole moment of the corresponding homogeneous (without 
solutes) polarized solvent and f2 is the volume of the solute; P = (47r) _1 (l — e^ 1 )Eo 
is the polarization of the homogeneous solvent. The second summand in equation (JSJ) 
represents the dipole moment cut from the liquid by inserting Nq voids. Finally, the 
last term is an additional polarization induced in the surrounding liquid by the surface 
charge density a p. 

The value of the interface solute dipole Mq 1 * will depend on the specifics of the 
solute-solvent interactions and the local polarization of the solvent created by these 
interactions. While it is a complex function of the entire mosaic of pairwise solute- 
solvent interactions for a realistic solute, an estimate of this parameter can be obtained 
from dielectric theories for a spherical void in a dielectric. The interface dipole reads in 
this case [15] 

M^ = -3fi P/(2e s + l), (6) 
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where the subscript "M" specifies Maxwell's electrostatics of a dividing surface not 
affected by local solute-solvent interactions. In order to quantify deviations from this 
generic result, one can introduce the ratio 

a = M'^/Mf. (7) 

The dipole moment of the mixture is related to the mixture dielectric constant e m i x 
as M mix /V = (47r) _1 (l — e~^ lx )E . One then obtains for the dielectric constant of the 
mixture 

e. 



-=- = l + ^(e,-l) l-2a-^-— +R x (vo), (8) 

where t/q = N Q /V is the volume fraction of the solutes in the sample with the overall 
volume V. We have put an extra term i?i(r/o) i n th e above equation to indicate terms 
non-linear in the volume fraction that appear in the dielectric constant when mutual 
polarization of the interfacial dipoles is taken into account [22] . Similar non-linear terms 
appear in the response of a mixture of water with dipolar solutes discussed below. There 
is presently no consistent formalism to include these effects and we neglect them at the 
current stage of the theory development recognizing that the theory might run into 
conflict with the data collected for concentrated solutions. 

If the Maxwell result for a void in a dielectric holds, a = 1 and the dielectric 
constant of the mixture becomes 

£ s -, 3(e s — 1) 

— = 1+ "»i-rr- (9) 



-mix 



Equation (jSJ), with Ri(rj ) omitted, and (|5J) describe the dielectric constant of an ideal 
mixture of non-polar solutes and a polar solvent. Equation (Q also reduces to the 
standard result of the Maxwell- Wagner theory in the limit of low volume fraction of the 
solutes J2J [3] • One can also account for the electronic polarizability of the protein not 
mentioned so far. If the refractive index n p can be assigned to the protein, one needs 
only to realize that the boundary conditions of the dielectric theories are sensitive to the 
ratio of the two dielectric constants at the dividing surface, e s /n^. Equation (jHJ) then 
extends to 

= l + Vo 9 I 2 • ( 10 ) 

Equation (JSJ) can be alternatively written in terms of the cavity field E c inside a 
spherical void in a uniformly polarized liquid. The electric field inside the cavity is 
proportional to the external field, with the susceptibility Xc = E c / Eq. In terms of this 
susceptibility, (jSJ) becomes [18] 

-^- = 1 + 3t/o [Xce s - 1] . (11) 

The standard prescription of Maxwell's theory of dielectrics predicts [HJ [23] 
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The connection between the susceptibility Xc and the parameter a (equation (j7j)) 
that is required to obtain equation (ITT]) from equation (|Sj) is derived from the following 
arguments. The polarization P(r) in the solvent, induced by the Maxwell field given by 
equation ([3]), creates a non- vanishing electric field inside the solute that is given by the 
equation 

E c = E + / T(r) • P(r)dr. (13) 



Upon substitution of equation (J3J) into this relation, one arrives at the connection 
between Xc and a 

2(e s - 1) 
3e sXc = e s + 2-a ^ + / - (14) 

Combining equations (J5J) and (113]) . one arrives at equation ( ITT]) . 

#.#. Polar Solutes in a Polar Solvent 

When a solute carries dipole moment mo, it aligns along the external field such that the 
average dipole {m )E in a weak external field is given by linear susceptibility [23] xo 

(m ) E = xo^oEo, (15) 

where (■ ■ -)e denotes an ensemble average in the presence of the external field and 

Xo = Xoo + Xos = ((3/3n )(5m ■ 5M mix ). (16) 

In this equation, 5m = m — (m } and 5M mix = M mix — (M mix ) are the deviations 
of the solute dipole and the dipole of the sample M m i x from their average values and 
(3 = 1/ik-QT) is the inverse temperature. 

The solute susceptibility in equation (Tl6]) is split into the self, Xoo> an d solute- 
solvent, Xos, parts. The former is given by the variance of a single solute dipole 

Xoo = (/3/3fi )((5m ) 2 ). (17) 

Correspondingly, the cross susceptibility is the correlation of a single solute dipole with 
the dipole moment 5M S of the entire solvent in the sample [21] 

Xos = ((3/3n )(5m -5M s }. (18) 

Equation (1T7]) neglects correlations between dipole moments of the solutes in the solution 
represented by the corresponding Kirkwood factor. Since the latter describes short-range 
correlations, of the length-scale of the molecular diameter [9] , they can be safely omitted 
in the type of theory developed here. 

Both standard arguments of the dielectric theories [H] and microscopic derivation 
[T2] suggest a simple connection between the solute dipolar susceptibility xo an d the 
self susceptibility % o 

Xo = XcXoo- (19) 

This relation implies that the account of the solute-solvent cross-correlations entering 
susceptibility xos amounts to introducing the cavity field acting on the average solute 
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dipoles, which also defines the torque acting on a selected dipole in the Onsager theory 
of dipolar liquids (directing field) [T3] . 

Adding the dipolar polarization of the solutes to equation (ITT]) for the dielectric 
constant of the liquid with spherical voids, one arrives at the dielectric constant of the 
solution 

= 1 -3r/o + 3?7oe s x c (l - j/o) , (20) 



where y = (47r/3)xoo- This equation clearly reduces to (ITT]) in the limit of non-polar 
solutes when y — > 0. 

2.3. Frequency- Dependent Response 

The static arguments presented in the previous sections can be extended to the frequency 
domain of main interest to broad-band dielectric spectroscopy. The dielectric constants 
of both the solvent and the mixture become frequency- dependent functions, e s (u) and 
e m i x (u;). The dipolar susceptibility of an isolated solute transforms into a linear response 
function, instead of a static correlator of equation (IT7|) . The relevant formalism is 
well developed and the result is the following response function of the solute dipolar 
fluctuations I2H 1251 



Xool^J = xoo 



1 + iu S 00 (u) 



(21) 



Here, Sqq(oj) is the Laplace- Fourier transform of the normalized time correlation function 
of the solute dipole m (t) 

S 00 (t) = [((Jmo) 2 )]" 1 (Sm (t) • <Jmo(0)>. (22) 

This function was fitted to multi-exponential decay when applied to the analysis of the 
MD simulation data presented below 

S m {t) = Y,^- t/T \ ^> = l, (23) 

i i 

where r» are the relaxation times and Ai are the relative weights of the relaxation 
components. From this equation, one gets the frequency-dependent function yo(u) 

y {u) = y Q Y jT ^— (24) 

i 

The frequency-dependent dielectric constant of the solution becomes 

^T = 1 - 377o + ^ Q e s {uj)xM (1 - 2/oM) . (25) 

^mix(,Wj 

Our arguments so far have not included any approximations except neglecting 
mutual polarization of solutes at their high concentration and the short-range 
correlations of solute dipoles entering the Kirkwood factor of the solutes. However, 
equations ( 12 Op and (123]) contain an unknown cavity-field susceptibility Xc(w). The 
Maxwell's result for this function refers to a free surface separating a dielectric from 
a void. It is a priory not obvious that this function can describe the complex and 



Protein dielectrics 



0.25 

0.2 

-o 0-15 

0.05 



-- HOmg/mL 


.-_ 


28 mg/mL 


■ /' 


57 mg/mL 


//" 


— Ill mg/mL 


/ 


water/buffer 


/ 


- 


/ 








i^ 


i ,i 


i , .... ...i 



0.01 



0.1 



1 10 

v, GHz 



100 



Figure 1: Real part of the cavity-field susceptibility Xc(w) extracted from experimental 
dielectric measurements according to (I25p . The results combine the broad-band 
dielectric measurements from [19] with high-frequency data from [20]. The dotted line 
indicates y^ from equation f[T2l for pure water (at lower frequencies) and the buffer (at 
higher frequencies). The gap between the two sets of curves represents the frequency 
window between the measurements. 



heterogeneous protein-water interface involving both weak protein-water interactions at 
hydrophobic patches and strong binding to charged surface residues. However, one can 
use the experimental input for the dielectric constants of the mixture and pure water in 
equation (1251) to extract the cavity-field susceptibility x c (w). 

Figure [T] shows the real part x' c ( u ) extracted from equation (l25l) using frequency- 
dependent dielectric constants of lysozyme solutions from broad-band dielectric 
spectroscopy below 50 GHz [19] and from separate measurements in the frequency range 
70-700 GHz [20]. The dotted line shows Re[xf] from equation (ED; the break in the 
curve signals the transition from water to buffer used at higher frequencies in [20] . The 
cavity-field susceptibility follows very closely the Maxwell prediction in the range of 
frequencies 10-200 GHz, but then deviates downward outside this range. The behavior 
at low frequencies is particularly noteworthy. 

It turns out that the dipole moment induced at the protein by an external field 
is over-screened [26] by the hydration layers, and perhaps the ionic atmosphere, to 
nearly zero. In fact, \ c is below zero at v < 1 GHz, implying a die-electric response, 
i.e. repulsion of the protein dipole from a region of a stronger electric field. This 
phenomenon, known as negative dielectrophoresis, is well-documented for hydrated 
nanoparticles [27], but has not been broadly observed for proteins. Our recent extensive 
simulations of ubiquitin [12], which is neutral at pH= 7.0, have indicated exactly this 
scenario: a negative xos, larger in magnitude than the positive x o, thus resulting in 
a slightly negative Xo i n equation (TlBl) . However, this result has not been detected by 
simulations of charged proteins, including lysozyme, probably due to the neglect of the 
ionic atmosphere in the analysis. 

Figure [1] suggests that dielectric models of the cavity-field susceptibility do not 
provide an adequate description in the entire range of frequencies of interest to broad- 
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band spectroscopy. However, the modeling can proceed along separate routes since the 
expulsion of polar water from the solute core is significant only at high frequencies, 
while the polar response of the protein dipole, described by yo(u), dominates at low 
frequencies. One therefore can keep the Maxwell result for x c (cj) for the former 
component, as realized in equation (Q. Since there is currently no model allowing to 
describe the overscreening observed at low frequencies, we have resorted to an empirical 
approximation. Replacing Xc^sVo in eqs (1201 and (j25|) with x^Vo accomplishes most of 
what is seen to occur in figure Q] and allows us to arrive at a compact relation for the 
dielectric constant of the solution 

1 + - f? - [e» - 1 - 3y H] ■ (26) 



e mix (w) 2e s (uj) 

2.4- Dielectric instability 

Equation ( |20l) predicts a point of dielectric instability at which the assumption of a 
uniform solution of weakly interacting protein dipoles breaks down. The instability 
is toward clustering of dipoles and is associated with the divergence of the dielectric 
constant e m ; x . It is reached at the critical volume fraction 

377 c =[H-e s Xc(2/o-l)] _1 . (27) 

If the Maxwell form of the cavity-field susceptibility is used in this equation, the critical 
point f] c = 0.01 (t/o — 16) corresponds to the concentration of 8 mg/mL for lysozyme 
in solution. Lysozyme solutions are stable in this range of concentrations and this 
estimate is clearly too low. On the contrary, the overscreening scenario shown in figure 
[1] makes r] c negative, thus removing the instability altogether. While other forms of 
aggregation are still possible [281 I2H], it might be quite possible that overscreening of 
the protein dipole eliminates instability toward dipolar clustering (such as formation 
of dipolar chains) and lowers the sensitivity of proteins in solutions to inhomogeneous 
electric fields always present in vivo. 

3. Application to Experiment: Lysozyme Solution 

Dielectric measurements of solutions typically provide the real and imaginary parts of the 
dielectric constant as functions of frequency and solution composition [301 EB Q21 E21 [33] . 
The existence of these two coordinates, frequency and solute concentration, allows one to 
learn about the specific pattern of interfacial polarization realized for a given solute and 
the dynamics of processes contributing to the relaxation of the sample dipole moment. 
We start with the analysis of the concentration dependence at a given frequency, followed 
with the analysis of the frequency dependence at a fixed concentration. 

3.1. Decrement of the water Debye peak 

Independently of the details of the dynamics of a protein itself and its coupling to 
the interfacial waters, the time-scales of these motions are significantly lower than the 
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Figure 2: Decrement of the water dielectric constant in the solution of lysozyme in 
water, Ae(u)/e s (u) = e mix (u)/e s (uj) — 1, as a function of the protein concentration C . 
The points are the experimental data at the frequency of the water Debye peak ud ~ 18 
GHz (circles) [19], and at v = 72 GHz (diamonds) [20]. The solid and dash-dotted 
lines refer to the calculations using equations Qj and ( TTOl) in the order of increasing 
frequency. The dashed line is the calculation incorporating the dynamics of the protein 
dipole according to equation ( 121)]) with yo(u>£>) calculated from MD simulations [T2]. The 
lysozyme molecular volume of f2 = 29.8 nm 3 is used to convert from the volume fraction 
to the solution concentration. The dotted line connects the experimental points. 



characteristic time of dielectric relaxation of water. The global motions of the solute 
are dynamically frozen at the frequency of the water Debye peak (z/q ~ 18 GHz). 
This implies that yo{uiD) can be dropped from equation ( )26|) . One then arrives at the 
dielectric constant of the mixture of polar water and effectively non-polar solutes (eqs 
OH]) and OH])). Any sufficiently high frequency can in principle be taken for this analysis. 
The decrement of the Debye peak of water in the solution vs the solute concentration 
is often reported [8] and can be used, in the framework of the present theory, as a 
convenient source of data to extract the information about the parameters a and \ c . 

Our formalism is applied to recent measurements of dielectric spectra of lysozyme 
solutions [19] [20]. Figure [2] shows the dependence of the decrement in the amplitude 
of the water Debye peak w D in the solution Ae(cjD) = emix(^D) ~ £ s (wd) vs the protein 
concentration. Circles show the experimental data from [19] , while the solid and dashed 
lines refer to equations ([9]) and (126]) . respectively. For the latter, t/oO-^d) calculated from 
MD simulations [12], and discussed below for the analysis at lower frequencies, was used. 
Clearly, the protein permanent dipole can be safely neglected. The transformation from 
the solution concentration reported experimentally to the volume fraction required by 
((9]) and ( |25l) was performed by using the volume of lysozyme Qq = 29.8 nm 3 . The latter 
was calculated from the crystallographic structure of the protein (3FE0, PBD database) 
by using the algorithm developed by Till and Ullmann [34] . 

In accord with the results shown in figure tU the cavity-field susceptibility is well 
described by the Maxwell form (equation ([12]) ) at the frequency of the water Debye 
peak, and the agreement between theory and experiment is excellent. It becomes less 
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Figure 3: Real (a) and imaginary (b) parts of the dielectric constant of the lysozyme 
solution measured experimentally [TJI| (points) and calculated theoretically (lines). The 
dotted lines show the real and imaginary parts of the dielectric spectrum of water from 



satisfactory at a higher frequency of 72 GHz [20] , also shown in figure [2j The refractive 
index of the protein starts to affect the result at this high frequency and n p = 1.7 from 
[35] was adopted in the calculations using equation (flOl) . The theoretical slope with 
increasing protein concentration is higher than experimentally reported and is likely 
related to deviations from the Maxwell form of the cavity-field susceptibility seen in 
figure [Q 



3.2. Dielectric spectra of solutions 

The results for the dispersion and loss spectra of lysozyme solutions are shown in 
figure [31 Experimental data from [36] were used for e s (u) and Molecular Dynamics 
(MD) simulations of a single lysozyme protein hydrated in a simulation box of TIP3P 
waters [12] were used to produce yo(oo) in ff2^1) . The relaxation parameters in (T24"]) are: 
Ai = {0.13,0.06,0.81}, n = {0.037,0.295, 14.6} ns, y = 16.3. The dominant relaxation 
component of the solute dipole, with the relaxation time of 14.6 ns, can be assigned 
to protein tumbling. The relaxation time of 9.1 ns was reported for this relaxation 
component from the analysis of proton NMR at low resonance frequencies |37j . 

The usefulness of MD simulations is somewhat limited for the sake of comparison 
with experiment since the charge distribution in the protein studied by simulations might 
not entirely fit the experimental conditions. The standard force-field prescriptions for 
protonating/deprotonating the surface residues of lysozyme at pH = 7.0 produce the 
overall protein charge of +7, while the charge of +10 is reported at pH= 5.5 in the 
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experimental study [TJ5]. Overall, permanent dipole moments of proteins arise from 
slight deviations from highly symmetric distribution of charge minimizing the total 
dipole moment [Ml EH]- Shifts of pK a values of surface residues due to local electrostatic 
environment [3D], ion association, and pH can therefore alter the dipole moment. 

Despite remaining uncertainties regarding the magnitude of the protein dipole when 
experimental conditions are concerned, the dipole (m ) = 223 D from MD results in 
a fair agreement between theoretical and experimental dispersion curves e' mix {ui) a ^ the 
lower concentration of the protein, 28 mg/mL (figure Ek)- The theory misses some of the 
static dielectric constant at the higher concentration, c = 110 mg/mL, but the difference 
actually comes from the missing increment at intermediate frequencies associated with 
5-dispersion. This part of the spectrum is also missing in the loss spectrum (figure [3b). 
This outcome is expected since specific protein-water binding contributing to this signal 
, I3D ED] has not been incorporated into the model. 



4. Summary 

Broad-band dielectric spectroscopy is a widely used tool to interrogate the dynamics 
of complex systems, including protein solutions. The interest in the field in the recent 
years has been to extract the polarization properties and dynamics of protein hydration 
layers from frequency- dependent spectra. The standard approach to the problem is to 
fit the dispersion and loss spectra to a sum of Debye or stretch-exponential functions, 
assuming that each component represents a separate relaxation process in a complex 
environment. The obvious limitation of this approach is the non-additivity of interfacial 
polarization, well recognized by classical theories of dielectric mixtures [2], El 13]- While 
these classical theories were developed for mixtures of dielectric materials, when each 
component can be represented by a macroscopic dielectric body, their application to 
hydrated proteins is clearly limited. At the same time, standard theories of dipolar 
liquids pUJ E] are not of much use either since they do not recognize the existence of 
an extended polarizable interface, which is in fact the central concept of the dielectric 
theories of mixtures. The present theoretical development aims to fill the void existing 
in each approach by recognizing both the molecular nature of the protein dipole and a 
quasi-macroscopic subensemble of interfacial waters producing interfacial polarization. 
The theory thus aims to study if the standard rules established for cavities carved 
in dielectrics, and also applied to calculate the local field acting on molecular dipoles 
[T3] . can be applied to hydrated proteins. Equations ( |25j) is central to this analysis 
since it allows us to extract the cavity-field susceptibility Xc(^) directly from the 
frequency- dependent dielectric constants of the protein solution and pure water. The 
remarkable result of this analysis is that at u < 1 GHz the susceptibility Xc{w) is below 
~ 0.02 predicted by the Maxwell equation f fl2|) and is in the negative territory, down to 
~ — 10 -3 . Therefore, the standard prescription derived for dielectric cavities (equations 
(JTJ) and (TT2l ) cannot be used in successful theories of dielectric response of protein 
solutions. 
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Granted, the cavity susceptibility extracted from experimental measurements might 
reflect the combined response of the dielectric interface and the ionic atmosphere. 
However, as a cumulative signature of the protein-water interface, it dramatically 
downscales the permanent dipole sensed by the dielectric experiment compared to its 
value calculated from atomic charges. Its low value can also help to explain the puzzling 
ability of proteins to stay in solution in vivo, despite significant electric field gradients 
that should pull a paraelectric particle to stick to, for instance, the bilipid membrane. 
The die-electric response suggested by the present analysis of experimental data, and 
our previous simulations (12], might be an answer to this puzzle since a die-electric 
solute repels from a charged interface creating the field gradient. It also eliminates the 
dielectric instability toward clustering of the solute dipoles predicted by (|20|) and (|2Tj) 
when the Maxwell form of the cavity-field susceptibility is used there. 
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